Integral relations for the surface transfer coefficients. 
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In this paper we derive relations between the local resistivities inside the interfacial region and 
the overall resistances of the surface between two phases for a mixture. These resistivities are the 
coefficients in the force-flux relations for the stationary heat and mass transfer through the inter- 
face. We have shown that interfacial resistances depend among other things on the enthalpy profile 
' across the interface. Since this variation is substantial (the enthalpy of evaporation is one of the 

main differences between liquid and vapor phases) the interfacial resistivities are also substantial. 
Particularly, surface put up much more resistance to the heat and mass transfer then the homoge- 
neous phase. This is the case not only for the pure heat conduction and diffusion but also for the 
cross effects like thermal diffusion. 

I. INTRODUCTION 

(N ' 

During evaporation and condensation heat and mass is transported through the interface. The common description 
of these phenomena uses certain assumptions which are debatable. For instance, one usually assumes equilibrium 
conditions at the interface [l| and neglects the coupling effects between the heat and mass transfer Q. The first 
one is a zero approximation and it is desirable to extend it to a more accurate theory. This is done in [2j for 
instance. Neglecting the coupling coefficients was shown to violate the second law of thermodynamics for a one- 
component system [3| . The coupling is important since the corresponding transfer coefficients depend on the enthalpy 
of vaporization. The significance of this quantity in this context is due to the large difference between the liquid and 
i-^j , the vapor values of the enthalpy. 

The situation becomes even more complicated when one consider evaporation in mixtures. These processes hap- 
pen in such industrial applications as distillation and therefore the precise description is important. Depending on 
. conditions one can get the mass fluxes of components in the same or in the opposite directions. 

There has been done a number of studies of the interfacial transport for one-c omp onent systems: experiments 
0, H, H, 0] , molecular dynamic simulations || [t| [T(| EH, E2 , kinetic theory [f3|, EH, EH, EH and square gradient 
continuous description [l?], EH- All these works use different approaches, which allows one to investigate different 
aspects. Mainly one-component systems have been studied. One of the points of interest is the dependence of the 
overall interfacial resistances on the continuous profiles. Once we have a description which relates the resistance of the 
Gibbs surface [I9j to the continuous profiles of, in particular, local resistivities, we can describe the surface separately 
from the adjacent bulk phases which is closely related to hypothesis of local equilibrium of a surface 0]. 

In this paper we extend the analysis done in [l8j for one-component system. Here we do the analysis for mixtures. 
Using the square gradient theory for the non-equilibrium interface developed earlier [20| we derive the general relations 
for the case of the transport perpendicular to the surface. We show that, as in the case of a one-component system, one 
can obtain the interfacial resistances using the continuous profiles obtained for equilibrium. This simplifies analysis a 
lot since one does not need to consider a non-equilibrium solution in order to obtain these resistivities. This is in fact 
a requirement, which the interfacial resistances must satisfy: as being defined within the linear force-flux relations 
jjj ■ they should depend only on unperturbed quantities, i.e. equilibrium ones. 

The evaporation and condensation often take place not only through planar interfaces, but also through curved 
ones, like the evaporation into a bubble. We do not restrict ourselves to the planar interfaces and give the expressions 
for the interfacial resistivities for curved surfaces. 

In [2l| we have obtained the overall interfacial resistivities R 9 and Ri using three different methods: an experiment- 
like procedure, a perturbation cell method and kinetic theory (only R 9 ). Those methods were found to be in a good 
agreement. In this paper we compare them with the interfacial resistivities found using integral relations, which 
relates these quantities to the local resistivities directly. We show that the agreement is also good. 

In this paper we will focus on the heat and mass transfer through the interface of a two-phase mixture. We will not 
consider non-equilibrium perturbation along the surface. We will also assume the fluid to be non- viscous. Furthermore 
we will focus on the stationary non-equilibrium perturbation. In Sec. [UJ we give the expressions for the local and 
the excess entropy production found in [2(| and [2f|. In Sec. [HI] we derive the integral relations in general form. 
It is convenient to use measurable heat and mass fluxes and we therefore show how to translate general relations to 
the resistivities associated with the measurable fluxes in Sec. |IV| . In Sec. [V] we give the explicit expressions for a 
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binary mixture and apply the analysis to the particular mixture of cyclohexane and n-hexane. Concluding remarks 
are given in Sec. | Vlj . 

II. THE ENTROPY PRODUCTION. 
A. Stationary states. 

Consider the total energy flux J e and the mass fluxes = Pi Vj and J m = pv, where pi and p are the density of 
i-th component and the overall density respectively, while Vj and v = J27=i are * ne velocity of the i-th component 
and the barycentric velocity respectively. Furthermore & is the mass fraction of i-th component and n is the total 
number of components. In the stationary states these fluxes satisfy the relations 

V-J e = 0, V-J^=0, V-J m = (II.l) 

As we are interested in transport through the surface we will restrict analysis to solutions of the form J(a;i,X||) = 
(J(xi), 0, 0), where J is one of the above fluxes. Furthermore x\ is the normal coordinate to the surface and X|| = 
(#2; £3) are the tangential coordinates. In the case of normal transport all the quantities depend only on x\ but not 
on X|| . For these solutions 

V-J = V ± J=— i— A(p, 2 f, 3 j) = (II.2) 
f)if) 2 f)3 dx 1 

where t)j = i)i(xi,x^) are Lame coefficients for curvilinear coordinates. It follows from Eq. (|II.2[) that 

f) 2 (zi) i) 3 (xi) J{xt) = t) 2 (* S ) M*") J(x s ) = |j« [)* J s (II.3) 

where x s is the chosen dividing surface and f)| and t)| are the Lame coefficients on that dividing surface. We will 
suppress the superscript s for the flux J as long as it does not lead to confusion. 

B. Local entropy production. 

Consider the local entropy production for the mixture interface found in [201 ] , written for the case of transport in 
the direction through the interface 

° s = J q - £ Ji Vxf = J q J2 J > ( IL4 ) 

i=l i=l 

where T is the temperature, J q is the heat flux and Ji = pi (v, — v) = — ^ J m is the diffusive mass flux, which 
satisfies the relation 



Ji - (II.5) 



i=i 



Furthermore ipi = pi — where /i^ is the chemical potential of the i-th component. The heat flux is related to the 
total energy flux and to the measurable heat flux J' q as 

n 

Jq = J e -J m (h + v 2 /2 - g-r) =Jg+J2 h -J ( IL6 ) 

i=l 

Here hi and h are the partial enthalpy of the i-th component and the specific enthalpy respectively, while g is the 
gravitational acceleration and v 2 /2 — g-r is the sum of specific kinetic and potential energies of the mixture. 

For the stationary transport through the interface the Gibbs-Duhem relation has the following form (see [2l[ for 
details): 



c dT 1 S^t dtli y 9an = Q (117) 

dx\ 1 dx\ dxi 

i—i 
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where an is an element of the thermodynamic pressure tensor cr a p. Together with the equation of motion, \7±<Jii 
pVi(v 2 /2 — g-r) = 0, it gives 



i—l 



1 



(ii. 



where fa = ^ + v 2 /2 - g-r and hi = h t + v 2 /2 — g-r. Substituting Eq. (|II.6[) and Jj into Eq. (|H.4|) and using Eq. (|II.8|) 
one can express the local entropy production (|II.4[) in terms of the total energy flux: 



n 

1 v— * /l 

a s = J e V_l— - V -L 



and in terms of the measurable heat flux: 

1 nl ( ib I 



i=l ^ 



(II.9) 



(11.10) 



where rji = hi — h n = hi — h n . Note that the difference between partial enthalpies with a tilde is equal to the difference 
without the tilde. In stationary states Eq. (|II.9p . Eq. (jll.lOp and Eq. pi.4[) are completely equivalent expressions. 



C. Excess entropy production. 

We define the excess 4>{x s ) of a density 4>{x s ) per unit of volume in curvilinear coordinates as 

1 



<t>(x s ) 



^2 ^3 JX 



dx 1 l) 1 i)2h3(t> ex (xi;x s ) 



where 



^(x^x 3 ) = (f>( Xl ) - (f> 9 ( Xl ) Q(x s - Xl ) - /(xi) 6(xi - x s ) 



(II-ll) 



(11.12) 



Here x s indicates the position of the chosen dividing surface while x 9,s and x e ' s are the boundaries of the interfacial 
region at the gas and liquid side respectively. These boundaries are chosen such that cf>(x 9 ' s ) = <j) 9 (x 9,s ) and tfi(x e ' s ) — 
cj) (x e,s ) with a certain accuracy. Superscripts I and g indicate the function <f> extrapolated from the liquid and gas to 
the surface region. We refer to [2l| for more extensive discussion. 

Taking the excess of the local entropy production given by Eq. (|II.9|) we obtain (see [2l[ for details) : 



<?s = Je. 



1 1 



~t ~g 

K_B_ 



(11.13) 



Here T e = T e (x s ) and T 9 = T 9 (x s ) are the temperatures extrapolated from the liquid and gas to the dividing surface. 
The analogous meaning have jlf and ji 9 . All the quantities, both the fluxes and the forces, arc evaluated at the 
dividing surface x s 



Using Eq. (|II.6|) we obtain 



and 



J'"[--- 



7r 



i=l 



i=l 



u e U 9 \ - ( 1 



1 



1 



(11.14) 



(11.15) 



where JJ' * and J g '' s are the measurable heat fluxes on the liquid and gas side respectively. Again Eq. (|II.14[) , Eq. (|II.15P 
and Eq. (|II.13|) are completely equivalent in stationary states. 
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III. INTEGRAL RELATIONS. 



Consider the entropy production Eq. pi.9[) . each term in which has a form J V^0. Following the common procedure 
in non-equilibrium thermodynamics one can write the force-flux relations for those entropy productions. Since all the 
terms have the same form, it is sufficient to consider only one force-flux pair. The phenomenological relation for that 
pair then reads 

V ± <f,(x 1 ) = ±^&=r(x 1 )J(x 1 ) (III.l) 
t)i axi 

for the local entropy production. The corresponding term in the excess entropy production in Eq. (|II.13p has a form 
J r ((fr — cf> 9 ). The phenomenological relation relation for this term reads 

^(x s )-(j)3{x s ) = R{x s ) J(x s ) (111.2) 

For the general case <f> and J must be replaced by a set as well as R and r by the corresponding matrix. 
Let us introduce excess operators £„ and € r which we will apply to a quantity <fi. Let 

x e ' s 

(£„{>} ee / tfciFji^rja;') (111.3) 



and 

x l ' s 

e r {0} ee fjf, (j| / d Xl -^- <p ex (r; x s ) (III.4) 

Jx 9 ' U2U3 

where <f> ex is defined by Eq. (|II.12[) . In cartesian coordinates the excess operators £„, € r and the excess ~, are given 
be the same expression when </> is a density per unit of volume. One should not confuse them however, as these 
operators have different meanings. The excess </> may be applied only to a volume density <fi and means the surface 
density. In contrast, neither £ n {0} nor (£ r {0} need to be applied to a volume density. In Eq. (|II.9j) <f> can be the 
inverse temperature or a chemical potential divided by the temperature. 
Applying <£„ operator to the both sides of Eq. (IIII.1|) one can show that 

<j> l (x s ) - <j) 9 {x s ) = € r {r} J(x s ) (HI.5) 

Comparing Eq. (|III.5|) with Eq. (|III.2|) we conclude that 

R = e r {r} (III.6) 

This is the general form of the integral relation for the resistivities r and R. Eq. pil.5[) together with Eq. pil.6[) are 
the most important and fundamental results of the paper. The generalization to a set of <j) and J is straightforward 
and will be applied in the rest of the paper. 

We then proceed to the explicit expressions for the phenological equations. The local entropy production (|II.9p 
produces the following force-flux relations 



T 

and the excess entropy production (|II.13[) produces 

1 1 



1 " 
V±- =r e qq J e -^2r^J Si 

(HI.7) 



i=i 



(Uli 



~l ~g 



i=l 



where J e and as well as T and fij are evaluated at the dividing surface x s . The off-diagonal coefficients of both 
sets satisfy the Onsager reciprocal relations. 
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Comparing Eq. (|III.7[) with Eq. (|III.8[) and using Eq. (|III.5[) together with Eq. (|III.6[) we may conclude that 

R e — g- [ r e \ 

1L qq *~r X' qq! 

R e qi =£r{r e qi } = <Er{rt q }=R! q (III.9) 
R]i =<£r{r? i } = <£ r {rf j } = Rt j 

Eq. (|III.9|) represents integral relations for the resistivity coefficients associated with the total energy flux. 



IV. MEASURABLE HEAT FLUXES. 



It is convenient to obtain the integral relations for the interface resistivities associated with the measurable heat 
fluxes rather then the total energy flux. We obtain the relations between resistivities, both local and for the whole 
surface. 

First a note regarding the dependence of the resistivities on the reference state of for instance the enthalpy. Both 
R e and r e coefficients depend on the reference state, as they are the resistivities associated with the absolute fluxes. 
The resistivities R' 9 and r', which will be defined below, are associated with the measurable fluxes and therefore 
independent of the reference state. 



A. The whole surface. 



In order to obtain the resistivities for the whole surface we consider Eq. (|II.14[) for the entropy production in terms 
of the measurable heat flux on the gas side of the surface. The analysis for the liquid side is completely equivalent to 
the one done for the gas side. 

The entropy production (|II.14jl produces the following phenomenological equations 

1 ™ 

— - R' 9 I'- 9 - R' 9 h 

'n 1 (IV.l) 







T' 


T3 


(H 




<(-~ 


—) 




T9 1 




T9 j 



The resistivities from Eq. pV.l[) and Eq. (|III.8|) are related, using Eq. ()II.6[) by 



R e =R' 9 

pe p ' 9 _i_ u9 p / g 

^qi — n-qi T 'H ^qq (IV.2) 

R% = Rj? + V Rj! + h 9 3 R' q ? + hf h 9 R> 9 



Eq. (jIV.ip and Eq. (|III.8[) are linear relations between the forces and the fluxes. In Eq. pV.2[) we should therefore 
use coexistence values of the enthalpies. 

Inverting Eq. (|IV.2[) and using the equilibrium enthalpies we obtain for the resistivities associated with the mea- 
surable heat flux 

R' 9 = R e 
L qq qq 

p ' 9 pe u9 pe 

n qi — ^qi 'H, eq ^qq (IV. 3) 

R' 9 = R e -h 9 R e -h 9 R e +Tl 9 Tl 9 R e 

J' 1 J l h eq JQ 3, eq L ql ' i, eq j, eq qq 



where hi tEq — hi, eq (x s ) — g-r 5 and all the quantities are evaluated at the chosen dividing surface. 
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B. Local resistivities. 



In order to write the linear laws for the local forces and fluxes we use Eq. (|II.10[) and Eq. (|II.4|) . The entropy 
production (|II.10p gives 



n-l 



(IV.4) 



(IV.5) 



i=l 

Ipi f 

V± y ~ ? = r Ji J 9 " 12 r 'n J * 

i=l 

while the entropy production (|II.4jl gives 

^ n—1 
Vj_ — = T qq J q — ^ ] Tqi Ji 

i=l 

Vj_ — = /"jq Jq — fjj Jj 

1=1 

The local resistivities r and r ' are related as follows 

r qq = r 99 

r qi = r ql -T li r q q (IV.6) 

= r j'i — 7 7i r j9 — % r 9» "T" ^i^J r qq 

Again, we should use the equilibrium profiles for r?j. 

We now relate the local resistivities associated with the measurable heat flux to the local resistivities associated 
with the total energy flux, similarly 1 to how it was done for the overall surface resistivities. Comparing Eq. (|IV.4|) 
with Eq. piljp we obtain 



r e _ 


r' 


99 


gq 


r e - = 

qi 


r' ~h 


r e = 

en 


r' ~h 

1 qq' l ri 



k=l 
n—1 



k=l 
n-l 



■qqhjh, - V Zk(rk e h t + r' k h 3 ) + hi + r^hj+ 

k~i (IV.7) 



n—1 n-l n—1 



T — r 
nn qq 



+ r ']i ~ 12 ^ r ki + r 3k) + 1212 r ki & &> 3, i = 1, n- 1 

k=l k=l 1=1 

n—1 n—1 n — 1 n — 1 

' m h n hi - ^2 &{ r Lhi + r' qk h n ) + r' q Ji n - ^ r'kii* + 12 12 r ki & i= 

fc=l fe=l fc=l Z=l 
n—1 n—1 n—1 

? >n "^E + r ,'fc) + 1212 T kl & 6 



fe=l fc=l i=l 



1 The details of this procedure are given in Appendix lAl . 
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and comparing Eq. (|IV.5[) with Eq. (|III.7|) we obtain 



11-1 



T qi ~ ' 11 



r 

.1' 



h — y]r q k£k +r qi , i = l,n-l 

k=l 
n-1 

k=l 

n-1 
fe=i 

n—1 n— 1 n— 1 

fc=i fe=i i=i 

n—1 n—1 n — 1 n — 1 

h 2 - h^2^ k {r ke + r qk ) + hr ql - ^ r ki £fc + X X rjw £ fc ^' * = n—1 

fe=l fc=l k=l 1=1 

n—1 n—1 n—1 

h 2 - h ^ £,k(r ke + r qk) + ^2^2 rkl ^ k £ 
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k=l k=l 1=1 

where h = h + v 2 /2 — g-r. 

Again, as in Eq. (|IV.2[) we should use the enthalpy and the mass fraction profiles in Eq. pV.7[) and Eq. pV.8[) in 
equilibrium. This leads to the relation hi = hi >eq — g-r and h = h eq — g-r. 

V. RESULTS FOR THE TWO COMPONENT MIXTURE. 

In this section we derive the integral relations for a binary mixture. The formula given in this section is not 
restricted to a particular binary mixture however. We do that for the sake of convenience as well as because we 
will apply them to a particular two component mixture. For simplicity of notation we will not write subscript eq for 
enthalpies, keeping in mind that all of them should be evaluated at equilibrium. 

The relations (|IV.7|) - (|IV.8|) between resistivities in case of binary mixture take the following form 



(V.l) 



r e = 


r' 






ii 


qq _ 






r ql ~ 


r qq hi 


+ r' qX 6 




r q2 ~ 


T-'qq h 2 


-r' q iti 




'11 — 


r' h 2 

1 qq <H 


+ 2r ql £ 2 hi + 


r[i& 


'21 — 


rqq h 2 


hi + r ql (£ 2 h 2 


-bhi) 


'22 — 


r' h 2 
' qq <H 


- 2r ql 6 h 2 + 





and 



qq „ 
r qq h + r q i & 



r„„ h-r q i£i 
\ 2 
i 2 

r 22 = r n h 2 - 2r ql £1 h + m £ 2 



qq ~ ' « X- » J _ (V2) 
rti = rqq h 2 + 2r ql 6 h + m & 



T 2 1 = r H h 2 + r q l{^ 2 - £l)/l - T n & £ 2 
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respectively. Using Eq. (|III.9p and Eq. PV.3P we therefore obtain 



and 



R' 9 = 
L qq 


£r {^ q } 








R ig _ 


£r {r' qq {h x - 








t>iq _ 

JX q2 — 










n'g _ 
K n — 


(f ir' (h, 

C r \ r qq \ni - 


ug\2 , 9 / t n rt,^ ug\ i 


r ll ?2/ 




n'9 _ 


<£r {r qq (hi- 




- K) - a (hi 




n'9 _ 
rl 22 — 


&r K g (/l2- 


- ftf ) 2 - 2r£ (fta - K) + 






R' 9 

x qq 


= (£ r {r qq } 








nig 


= £r{r qq (h 


-M)+r,i&} 






Kl 


= £ r {r qq (h 


- /if) - r q i £1} 






r>'g 
-"-11 


= £ r {r qq (h 


-/ l ?) 2 + 2r, 1 6 (fc - /if) + 






r' 9 

n 12 


= £r{r qq (h 


-h{){h-h a 2 )+r ql {& (h- 


- K) -ti(h- 


/«?))-m66} 


R' 9 

lx 22 


= <E r {r qq (h- 


-h 9 ) 2 -2r ql ^ (h-h 9 ) + 







(V.3) 



(V.4) 



Furthermore we used that hi — h 9 — hi — h 9 in both equations. Eq. (|V.3[) and Eq. (|V.4[) are the integral relations for 
the resistivities associated with the measurable heat flux for a binary mixture. They clearly do not depend on the 
reference chosen for the enthalpies. 

We further need the values for the local resistivities. In [12] we developed a scheme which uses r coefficients in 
calculation. We therefore use the following expressions 



r qq {x) 


— r g _ 
qq 




-r 9 q)qo(x) - 


1" a qq ( r qq 


1- r 9 qq ) qi(x) 






= r ql ~ 




-r ql )qo{x) - 




^ r ql)n( X ) 


(V.5) 


rn{x) 


= r ll ~ 




- r fi) Qo(x) - 


f- an (r^ - 


\-rn)qi(x) 





where qo(x) and qi(x) for each resistivity are modulatory curves for the resistivity profiles: qo(x) is a smooth arctan- 
like function which changes its value from to 1 within the range [x 9 ' s ; x e ' s ] and qi(x) is zero on the boundaries of the 
[x 9:S ; x £:S ] interval and has a peak proportional to the square gradient of the density inside this interval. Thus, the 
first two terms in each expression for the resistivity represents a smooth transitions from the gas bulk resistivity to 
the liquid bulk resistivity, while the third term represents a peak in the resistivity proportional to the square gradient 
of the density. 

The bulk values r l and r 9 are related to the measurable transport coefficients such as heat conductivity, the diffusion 
coefficient and the Soret coefficient. We refer for the details to p2l |. 

Consider a binary mixture of cyclohexane and n-hexane, which is the same as in p2l | and [2lj |. Furthermore we 
consider a planar interface between liquid and vapor. The mixture is in a box with gravity directed along axes x from 
left to right. The gas phase is therefore in the left part of the box and the liquid is in the right part. 

We compare the resistivities found in [2l| to the values obtained from Eq. (|V.4|) . The relative difference be- 
tween them is almost the same within the range of temperatures and chemical potential differences considered: 
T = {325, • • • , 335} and ip — {400, • • • , 1000}. In Table [I] we give the relative errors for gas and liquid side coeffi- 
cients. The data are obtained for a-amplitudes, values of which were found to fit best kinetic theory in [2l| . 
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TABLE I: Relative error in percent between the gas- and liquid- side co- 
efficients obtained by "perturbation cell" and "integral relations" meth- 





ods at T eq = 


330 and ij) e q 


= 700 for f3 


= 0.0002 and 


Otqq — 9, OL\q 


= 0, 




an = 3 . 












phase 


Rqq 


Ru 


R22 


R q i 


Rq2 


R12 


gas 


0.019090 


0.064642 


0.058851 


0.020649 


0.020680 


0.097096 


liquid 


0.019090 


0.006266 


0.000432 


0.036270 


0.034886 


6.233983 



The relative differences are not more then a few promille. It is larger only for R l l2 which is discussed below. 
We also do a consistency check. Consider Eq. (jIV.ljl for two component mixture, which has a form 

X q = Rqq Jq — R' ql Jfj — R' q2 J^ 2 

Xi = R{ g J q ' - R{l J^ ~ R{ 2 (V.6) 

X% = R' 2q Jq - i?2! % - R' 22 ^2 

The left hand side of each equation must be equal to the right hand side. The difference therefore reflects the error. 
We give the relative error between the left an the right hand side of Eq. (jV.6|) in percent in Table pi] . As a testing 
perturbation we used one of those used in the perturbation cell method. 

TABLE II: Relative error in percent between the left- and right- hand 
side of Eq. (|V.6|I for coefficients obtained by "perturbation cell" and 
"integral relations" methods at T eq = 330 and Ve« = 700 for (3 = 0.0002 
and a q q — 9, ai q = 0, an = 3 . 







Integral relations 




Perturbation cell 


phase 


Xq 


Xi X2 


Xq 


Xi X2 


gas 
liquid 


0.059489 
0.059489 


0.037918 0.296959 
0.172608 0.027275 


0.046965 
0.046851 


0.087411 0.867098 
0.216819 0.014248 



Again, the relative difference is not more then a few promille. Given that this is the case even for a few percent 
difference in one of the coefficients, we may conclude that the values of the forces are insensitive to the precise value 
of this resistivity coefficient. This also indicates that the value of this coefficient obtained in [2l[ has a 6% error. This 
does not necessarily affect, however, the accuracy of the integral relations. 



VI. DISCUSSION AND CONCLUSIONS. 



In this paper we have derived integral relations for the resistivities to the transport of heat and mass across 
the interface for mixtures. We have given relations between the local resistivity profiles and the overall interfacial 
resistivities. 

The integral relations make it possible to calculate the interfacial resistivities in a relatively simple way using only 
the equilibrium profiles of the system. This is important especially for mixtures, for which the computation of a 
non-equilibrium profiles is much more time consuming. 

The integral relations give an insight in the origin of the interfacial resistances. According to Eq. (|V.3[) and Eq. (|V.4[) 
the interfacial resistivities depend on the variation of the enthalpy across the interface. The transport coefficients 
depend on the equilibrium enthalpies which vary a lot through the interface. One can see from the above formulae, 
that the dependence on the enthalpy of evaporation (the difference between the enthalpies of the liquid and gas phases) 
is crucial not only for the diagonal diffusion coefficient, but also for the off-diagonal coefficients. This is aa important 
result since cross coefficients are usually neglected in the description of the interfacial phenomena. 

Another factor which affects the overall interfacial resistivity is the local resistivity profile. For instance, for the 
heat resistivity this is the only factor. It is noticeable that the interfacial resistivity depends on the whole profile of 
the local resistivity, not only on its bulk values. It is therefore crucial to have complete information about the local 
resistivity profiles. We have used sums of a function that smoothly connects the liquid and the vapor values and 
a peak proportional to the square gradient of the density. In principle one can use any model for this and further 
investigation are required. We have shown in 21| by comparison with the predictions from kinetic theory, that the 
local resistivities do have a peak in their profile and that the overall resistivities therefore depend on the amplitudes 
of these peaks. Within the current theory these amplitudes are adjustable parameters. 

The integral relations are in fact mathematical equalities. Given the local resistivities defined through the local 
force- flux relations, for instance r e from Eq. (|III.7[) . Eq. 8[) follows and one can consider Eq. (|V.3[) as a definition 
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of the overall interfacial resistivities R e used in Eq. (|III.8|) . It means that the force-flux relations (|IV.1[) for the whole 
surface follow from the local force-flux relations. One therefore does not need the excess entropy production Eq. (|II.14p 
to obtain Eq. (|TTTT8|) . 

This allows us to use them as a test for the accuracy of the numerical solution of the non-equilibrium case. Given 
that both local and overall linear laws are true independently, the different methods to obtain the overall resistivities 
give information about the accuracy of the method. The discussion below Table [I] and Table [IT] is based on this 
observation. 



APPENDIX A: LOCAL RESISTIVITIES. 



We need to relate the resistivities r from Eq. (|III.7j) to the resistivities p from Eq. (|IV.4j) . This is done by comparing 
the coefficients at the same fluxes in these equations. To do this we need to translate the set of fluxes used in Eq. (|IV.4p . 
{J', Ji, • • • , J n -i}j to the set of fluxes used in Eq. (|III.7|) . {J e , • • • , J^„}- This is done with the help of the relation 



J i - J& ~ & E J i> 

k=l 

n 

J q = J e - E h kJ^ k 



(A.l) 



k=l 



Substituting J' and Ji into the first of Eq. (liV.4|) we obtain 

n— 1 n— 1 



V ^ f = r 'qq Je " E J t> ( r ?<A + ~ E r 'qktk) ~ J U {j'qqK ~ E r 'qk£k) ( A ' 2 ) 



i—k i—k 

Comparing it with the first of Eq. T[) we obtain 



T — T 
qq qq 

n-1 



•qqht ~ E r 'qk & + T 'qii t = 1, n-1 

fc=l 
n-1 

Tqqhn ~ E T 'qk & 



(A.3) 



k=l 

which are the first 3 equations of Eq. (|IV.7[) . 

In order to obtain the remaining relations we consider the second of Eq. pil.7[) , which gives 

E & v y ■'■ E -E J c.E <l<ti 

3=1 j=l i=l j=l 

~ n 

V ± ^-V ± >^ =J e (r? q -r° e )-J2j( i (.r] i -r e ni ), j = 1^1 



(A.4) 



Furthermore we use Eq. pi.8[) which in case of the transport in the direction only perpendicular to the surface 
becomes 

f>(v±£-fciVxi)=0 (A.5) 



Together with the second of Eq. pV.4[) it gives 

n — n ^ 

E^ v ^ = E^* V 4 



1=1 i ' - 1 (A - 6) 
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Substituting Vj_(l/T) from Eq. (|A.2|) and J' and Ji from Eq. (|A. 1|) we obtain the left hand size of Eq. (|A.6[) expressed 
in terms of the fluxes J e and and the resistivities r'. Comparing the result with Eq. (|A.4[) we obtain the following 
equations sets 

E"=i rjqtj = r q q ELi & ^* 

(A.7a) 

r jij — r ne = r qq r lj + r j<j, J = 1>71— 1 



E"=l = ( r qq k n ~ ELl &) ELl & ^ 

(A.7b) 

r f„ - C = O^j + ^'g)/^ - ELl ( r g'fe % + r jfc) &> 3 = i, n- 1 



E"=l r jA = i^qq^i ~ ELl r g'fe £fc + r g'i) ELl & fe fc> i = 1, 71-1 

(A.7c) 

*ji ~ r ni = (rqqVj + r- g )hi - ELi 1 ( r qk Vj + r jk) & + ( r g'i »7j + r-i), 3 = 1, n-1, » = 1, n-1 

solving which we obtain the relations between the remaining resistivities 



r e = 

3Q 


n-1 




fc=l 
n-1 


r e = 

ne 






k=l 
n-1 


r e . = 

J 1 


rqghjhi - ^6=0 
fe=l 



fe=i fc=i fe=i z=i 

n — 1 7i—l n — 1 7i — 1 

•qqhjhn - ^ ik(r' k Ji n + r' qk hj) + r/,ft„ - r jk £k + X) 12 r kl €l, 3 = l,n-l 
fe=i fc=i fe=i i=i 

n — 1 7i—l n— In— 1 

•gg/i„/ii - X tk(rLhi + r 'qkK) + r' q jtin ~Y r ^^ + YYl r ki & 6. * = I,"-! 



fc=l fe=l fe=l (=1 

n — 1 n— 1 n— 1 

r'nn = r' qq h 2 n - h n £ U< e + r qk) + J2J2 r ki t k & 

k=l k=l 1=1 



As one can confirm the symmetry of the r '-matrix leads to the symmetry of the r e -matrix and vice versa. We 
therefore do not give the expressions for rf_, and rj n in Eq. (|IV.7|) . 

The relations (|IV.8[) between the r e - and r- resistivities are derived in the similar manner. 
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